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Abstract —This paper proposes a relay selection scheme that 
aims to improve the end-to-end symbol error rate (SER) perfor¬ 
mance of a two-way relay network (TWRN). The TWRN consists 
of two single-antenna sources and multiple relays employing 
decode-and-forward (DF) protocol. It is shown that the SER 
performance is determined by the minimum decision distance 
(DD) observed in the TWRN. However, the minimum DD is likely 
to be made arbitrarily small by channel fading. To tackle this 
problem, a phase rotation (PR) aided relay selection (RS) scheme 
is proposed to enlarge the minium DD, which in turn improves 
the SER performance. The proposed PR based scheme rotates 
the phases of the transmitted symbols of one source and of the 
selected relay according to the channel state information, aiming 
for increasing all DDs to be above a desired bound. The lower 
bound is further optimized by using a MaxMin-RS criterion 
associated with the channel gains. It is demonstrated that the PR 
aided MaxMin-RS approach achieves full diversity gain and an 
improved array gain. Furthermore, compared with the existing 
DF based schemes, the proposed scheme allows more flexible 
relay antenna configurations. 

Index Terms —Decode-and-forward, beamforming, relay selec¬ 
tion, network coding, MIMO, full diversity. 


I. Introduction 

T Wo-way relaying (TWR) is a promising technique to 
improve the coverage and connectivity of relay aided 
networks □-0I. In a typical TWR channel (TWRC) 0, 
0, two source nodes exchange information simultaneously 
with the aid of a relay node. Assuming the absence of a 
direct link between the two source nodes, communication 
takes place in two stages: the multiple access (MA) stage 
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and the broadcast (BC) stage. During the MA stage, both 
source nodes transmit their individual signals to the relay 
node simultaneously. Then, the relay node broadcasts the 
processed signals to the source nodes in the BC stage. It 
is also assumed that the TWR schemes are assisted by the 
network coding in analog or digital domain 0, Q. In fading 
wireless channels, multiple relay antennas can bring diversity 
gain to the TWR system 0. This paper focuses on achieving 
diversity in the network coding aided TWR systems composed 
of multiple relay antennas and two single-antenna sources. The 
performance metric considered is the symbol error rate (SER), 
based on which the achievable diversity gain is derived. 

In terms of the symbol error rate (SER) performance, 
diversity techniques have been studied extensively for TWR 
systems. These existing works generally obtain full diversity 
gain by processing signal(s) from either all relay antennas or 
a single one selected relay antenna. To be more specific, in 
03 , 03 , the received signals from all relay antennas are 
combined, which allows both diversity gain and array gain 
to be achieved. The signal combining operation employed in 
03 , 03 obstructs its application in the distributed scenario, 
where the multiple relay antennas are distributed among sev¬ 
eral single-antenna relays. This limitation is removed by the 
distributed-antenna space-time coding (DSTC) based scheme 
m and a variety of antenna selection (AS) based schemes 
(10|— (121. In particular, the DSTC is implemented in 0 
with the cooperation among all antennas. The performance 
achieved by the scheme of 0 is further enhanced by the 
AS based techniques fT0l - fT2l . All these AS based schemes 
09-03 use the MaxMin criterion, which maximizes the 
quality indicator of the worst link originating from the selected 
relay to the two source nodes. The MaxMin criterion is able 
to obtain full diversity gain by employing the amplify-and- 
forward (AF) relaying protocol to the selected antenna HOlK 
DD. The decode-and-forward (DF) protocol is considered 
by a opportunistic two-way relaying (O-TR) scheme (H, 
which also invokes MaxMin criterion and the resultant SER 
performance exhibits full diversity order, although the scheme 
is aided by perfect error-correction-code (ECC). In other 
words, the DF protocol in the O-TR scheme is carried out 
on the codeword by codeword basis. However, it is observed 
that in the absence of ECC, the MaxMin criterion aided DF 
protocol no longer guarantees full diversity gain. 

Against this background, in this paper we aim to tackle 
the problem of full-diversity guaranteed uncoded transmission, 
where the DF protocol is implemented on the symbol by 
symbol basis. We focus on the scenario with K (K > 1) 


2 


relays each equipped with an arbitrary number of antennas. 
Hence, the existing scenarios considered in |9l- fl5l are its 
special cases. Intuitively, a natural transmission approach is the 
straightforward extension of the previous AS based schemes, 
which always select a single antenna despite the distribution 
of the total relaying antennas mounted on K relay(s). By 
comparison, instead of simply extending these works, in this 
paper we propose a relay selection (RS) based scheme, which 
not only gleans full diversity gain from relay selection, but 
also improves array array power gain by combining signals 
from multiple antennas of the selected relay. Therefore, the 
proposed RS scheme enjoys the advantage of array array 
power gain over the existing AS based schemes. 

More specifically, the key enabler of the proposed RS 
scheme is a phase rotation (PR) strategy. In general, the 
PR strategy is executed symbol by symbol, which intends 
to remove the randomness of the phase metric of TWRC. 
This phase metric is termed effective angle (EA). In our 
previous work fRl . the randomness of EA is shown to 
degrade the achievable diversity gain in the MA stage. The 
proposed scheme in lfl4l is applicable to either the single 
multi-antenna relay or the selected single-antenna relay, and 
it employs a sub-optimal detector specific for the MPSK 
modulation. By contrast, in this paper we consider a more 
flexible antenna distribution and tackle the EA randomness 
problem under arbitrary modulation. Throughout this paper, 
the selected relay employs the optimal maximum likelihood 
(ML) detector, which is in principle more fundamental than 
the sub-optimal detector proposed in M. Furthermore, we 
point out that in the more general scenario considered, the 
randomness of EA impose an impact on the achievable per¬ 
formance in both the MA and BC stages. Correspondingly, as 
the derivatives of the PR strategy, the PR preprocessing and 
PR beamforming operations are conceived for the MA and 
BC stages, respectively. By removing the randomness of EA 
in both stages, the PR operations shape the distribution of the 
decision-distances (DDs), i.e., the Euclidean distances between 
all desired received signal points, as observed by receivers. It 
is demonstrated that the stochastic property of reshaped DDs 
allows the MaxMin-RS criterion to attain full diversity when 
the DF protocol is executed symbol by symbol. 

Notice that some existing approaches also process phase of 
symbols transmitted in one-way relay system which seems 
a little bit similar to the proposed PR strategy na-na. 
However, the full diversity gain achieved by these approaches 
is measured by outage capacity. This paper focuses on the 
influence of PR strategy on SER performance. All the diversity 
analysis provided by this paper are in the sense of SER perfor¬ 
mance rather than outage capacity. This is striking difference 
between our work and the other existing paper na-na. 
Generally, the diversity analysis method and its result derived 
from outage capacity are not equivalent to that from SER 
performance. For example, it is demonstrated that according 
to outage capacity, MaxMin-RS without any preprocessing is 
capable of obtaining full diversity gain fT5l . However, our 
previous work has shown SER performance does not exhibit 
any diversity gain by applying MaxMin-RS straightforwardly 
Da. This paper indeed provides insights into how PR strategy 


enhance the capability of MaxMin-RS on achieving diversity 
gains in SER performance. 

In particular, we provide theoretical analysis to confirm the 
impact of pre-canceling the randomness of the EA on achiev¬ 
able SER performance. The full diversity order property and 
the achievable SER performance of the proposed scheme are 
explicitly analyzed. Our analysis indicates that by exploiting 
the information embedded in the TWRC’s phase metric, the 
SER performance can be improved by using the MaxMin-RS 
approach in uncoded DF based relay networks. 

The rest of this paper is organized as follows. In Section II, 
the system model is described. In Section III, we present the 
proposed PR-MaxMin-RS scheme. In Section IV, we provide 
more details of the PR based approach, and derive the SER 
expression under certain system configurations. In Section V, 
numerical results are offered to confirm the advantages of the 
proposed scheme and to validate the theoretical analysis of the 
diversity order. Finally, we conclude the paper in Section VI. 

Notation : For a matrix, (-) H ,{-) T , and ||-|| represent the 
conjugate transpose, transpose, and Frobenius norm of the 
matrix, respectively. For a complex-valued variable, Sft(-), 
9 (•), | • |, and (•)* denote the real part, the imaginary part, the 
absolute value and the conjugate of the complex-valued vari¬ 
able, respectively. For a set, | • | represents the size of this set. 
CA/"(0,K) denotes a circularly symmetric Gaussian random 
vector with mean 0 and covariance matrix K. J\f (0,1) denotes 
a standard Gaussian distribution. Z- stands for the angle of a 
complex-valued number. 0 mxN denotes the M x TV matrix 
whose all entries are zero, and I mxN denotes the M x N unit 
matrix. E {•} represents the expectation operation with respect 
to its argument. Q (•) is the tail of the probability density 
function (PDF) of a standard Gaussian random variable, i.e., 

Q ( x ) = h Jo f ex P (-2 ot) de - 

II. SYSTEM MODEL 

A. Channel Model 

We consider the DF based TWRC shown in Fig.l, where 
two single-antenna source nodes S* (i = 1,2) exchange mes¬ 
sages with the aid of K half-duplex relay nodes denoted as 
{R/c |k = 1,..., K}. The relay node Rk is equipped with Lj , c 

K 

antennas, where Lk > 1, for k = 1,..., K. Let L = ^ Lk. 

k=l 

We assume that there is no direct link connecting the source 
nodes, while the links between the two source nodes and 
relays suffer small-scale fading. Specifically, the links from 
Si and S 2 to Rk are characterized by the x 1 channel 
vectors and g^, respectively. The elements hki and gki 
(l = 1, • • • , Lk ) of the vectors and are independent and 
identically distributed CAT(0,1). Phase offsets between carrier 
frequencies employed by the sources and Rk are absorbed into 
hki and g^i . It is also assumed that the channels are reciprocal 
(this can be true when using time-division duplexing (TDD)) 
and static during the two transmission stages for transmitting 
numerous data packets consecutively (this is true for block¬ 
fading channels). As a result, a small fraction of time used 
for the training symbols of CSI estimation and for feeding 
back the estimated CSI may essentially have little impact 
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Fig. 1. The TWRC system model considered: only one relay is selected and 
activated throughout the MA and BC stages. 


on the achievable transmission rate. At the beginning of the 
proposed scheme, one relay is selected to assist the source 
nodes in communication. The RS process can be implemented 
by using a distributed mechanism (T 9 ), which requires the 
relay node Rk to know only and g^, rather than the 
CSI of the entire network. Generally, the source nodes need 
extra feedback information from the relay to preprocess their 
information. This requirement can be removed under some 
specific configurations, which will be detailed in the sequel. 

We assume that Si and S2 share the same M-ary constel¬ 
lation alphabet set S, where \S\ = M. For the source S*, i = 
1 , 2 , the modulated symbol Si E S carries log M bits which are 
stacked in a 1 x log M vector Wi . The symbol Si is transmitted 
with power P s in the form of xi = y/PuiSi, where Ui denotes 
the preprocessing imposed on Si, and we have \ui\ = 1 . In 
addition, we define the mapping function M by which Wi 
is mapped to i.e., = Ml (Wi). Based on all, we define 

symbol-level XOR as si O s 2 = M (Ml ~ 1 (si) © Af -1 (s 2 )), 
where © denotes the bitwise XOR operation, and is the 
inverse mapping of Ml. Then, network coded symbol (NCS) 
is defined as snc — si © s 2 , whose alphabet is also S. 

B. DF Based Transmission with an Arbitrary Relay Rk 

We proceed to outline the transmission flow of the DF based 
TWRC. We consider the communication using an arbitrary 
relay Rk as an example. First of all, if preprocessing is 
employed by the sources, the relay Rk calculates and feeds 
back u\ and u 2 to Si and S2, respectively. The feedback is 
assumed to be perfect. Then, the bidirectional transmission 
will take place in the MA and BC stages. In the MA stage, 
the signals transmitted by the sources are assumed to arrive at 
Rk simultaneously. The signal vector received at Rk is given 
by 

yfc = h k xi + g k X2 + n, (1) 

where n is an Lk x 1 complex Gaussian noise vector obeying 
CM ( 0 Lk xi,cr 2 lL k xL k ). The ML multi-user detector (MUD) is 
employed to jointly decode both messages from y^, i.e., 

(S1.S2) 


arg min 

(?i ,fa)eSxS 


yfc - 


V Ps^-k'U 101 s/ ■P s§&;^'202 


2 

( 2 ) 


Then, Rk generates the NCS by using snc = si © £2, $nc © 
S. During the BC stage, Rk broadcasts the signal X# = 
VW^kSNC to the sources, where P r is the broadcast power, 
Wk is an Lk x 1 beamforming vector designed to exploit the 
spatial diversity offered by the multiple antennas of R&, and 

ll w fc|| = !• 

The signals received by Si and S2 may be expressed as 

Vbc,i = w fc sjvc + ni, 

VBC ,2 = \JFgk w kSNC + ri 2 , ( 3 ) 


respectively, where n\ and are independent Gaussian noises 
following CM (0, cr 2 ). In order to estimate snc , Si and S2 
employ the following ML detectors 


snc,i = arg mm 
es 

snc,2 = arg min 
<t>2Z.S 


Ubc, 


1 - 


Vbc, 2 


- CFgl w k<i>2 


( 4 ) 


respectively. For i = 1 , 2 , the source S3_* decides the desired 
information sent by according to Si = SNC,3-i 0 S3-1. 


C. The SER Bound When Using Rk 

For the relay Rk , the instantaneous overall SER serves 
as a relevant performance metric for the RS scheme dis¬ 
cussed later. Hence, let us analyze the instantaneous over- 
all SER, i.e., pf 2B 4 I £?=i Pr{S; ^ * |g fc ,h fc }. Accord¬ 
ing to EQl, we have pf 2£ < p^ A + \ EiLi pf f, 
where p^ A = Pr{sNc ^ snc |g*, h fc } and p^f = 
Pr jsjvc,i 7^ snc gA-- li/, j represent the SERs at the relay 
and the two sources, respectively. Furthermore, again rely¬ 
ing on ( 20 ), Pk 2E approaches + \ Y?i=i pf,f as the 
signal-to-noise ratio (SNR) increases. Thus, we use + 
\ Ei=i P k? approximate p E2E in the high-SNR scenario. 
In the sequel, we will analyze p E2E by first deriving the SER 
bounds for p^ A and pff. Then, the relevant performance 
metric will be extracted to enable the RS methods. 

We assume that when (0 i, 02) is transmitted the MUD ob¬ 
tains the estimate ^ 0' l5 4 >' 2 j of 0 i © 02 7^ 0 i 0 02 , 0 *, © S. 

Referring to snc — «i © $2, we know that the error event 
{s nc 7^ snc} comprises all possible cases. Based on this 
observation, pjf A can be given by 


MA 

Pfc 


— Pr i _ _ ©L _ ( 01 , 0 2 ) ^ ( 0 i, 02 ) |gfe,h fc l . 

( 01©02 7 ^ 0 'i© 02 ’ ' / V / I 


0Then, the union bound is used to approximate Wk ~ as 


^MA 


MA 

Pfc 


< 



1 

M2 




1 (01,02) -©■ 
transmitted, 


>i 5 0 2 ) denotes the event that when 
is obtained as its estimate. 


(01,02^ 


is 
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• Q (^\j^\\gkU2 (j>2 ~ <j>2^ + hfcMl (</>l - 4>'l) ||^ (6) 

(di,d 2 )G£) VV 7 W / 

(7) 

where p = |f is the SNR, D = 

|(dl,d 2 ) |rfi = <t>i- 0i,0 1 ©02 7 ^ 01 0 02, 0i,0i G 5 }, 
A(i/idi,i/ 2 d 2 ) — ||g/c^ 2 ^i + hfeiiid 2 || is the DD in 

the MA stage, and Alb denotes the lower bound of 
A (uidi, i/ 2 d 2 ),(di, d 2 ) G D. The first inequality given in 
© is based on the union bound whose tightness is well 
confirmed by J2T), while the tightness of the second upper 
bound in © depends on the design of A lb, which will be 
given later. 

On the other hand, the SER expressions in the BC stage 
conditioned on the instantaneous channel state are given by 

Pm = CxQ (ycmi |w^hfc| 2 ^ , 

P k,2 =CiQ (^C 2 pn |w£g fc | 2 ^ , (8) 


where p = |g, C\ and C 2 are constants depending on the 

modulation. Generally, and \j C 2 \w^g k \ 2 are 

DDs in BC stage. 

From © and ®, pf : 2E < + \ Ya=\ pf f is bounded 

[(!)]] by 



where a = + C\ is a modulation-specific constant. Since 

the accuracy of ® is ensured by the existing performance 
analysis j2TI . its substitution into p E2E < p^ A -\-^ Y^= i Pf,f 
does not impact the tightness of p E2E - Hence, the tightness 
of the upper bound given in © depends only on the design of 
A lb- We will propose a preprocessing scheme in Section III 
in order to find a desired A lb- Upon using the obtained A lb, 
the upper bound given by the right-hand side of © exhibits 
a reasonably small gap@, especially in high-SNR region, from 
the actual p E2E - As a beneficial result, the upper bound of 
p E2E in © may be used as a relay selection metric for 
achieving full diversity gain. 

More specifically, we propose to select the kth relay by 
minimizing the upper bound of p E2E which is given by ©. 
Then, we have 


k = arg max min 

k=l,...,K 


\2 i i 

a LB ri T, 

—, C2P w fc h 


C 2 P Wfcgfc 



This RS approach is not optimal. Because in each instant 
the criterion (ITOl) optimizes the upper bound of the SER 
performance, rather than the actual SER performance. The 
goal of this paper is to find a simple suboptimal RS approach 
capable of achieving full diversity gain. In the next section, we 
will design m, 112 and in order to ensure that the tightness 


2 Due to page limitations, numerical results illustrating this gap are not 
included in this paper. 


of the upper bound given in © does not impose a negative 
impact on the achievable diversity performance. The designs 
on ui, U 2 and are based on the PR strategy. As a result, 
the RS approach specified by (ITOl) is actually the MaxMin-RS 
with respect to channel gains. The PR strategy conceived not 
only assists the MaxMin-RS in achieving full diversity, but 
also allows us to obtain explicit analytical results. Notably, 
our analysis confirms that the proposed scheme obtains full 
diversity gain, and achieves better array power gain than the 
existing schemes. 

III. The PR AIDED MAXMIN-RS Scheme 
In this section, we first explain how the PR strategy shapes 
A (uidi, i/id 2 ), |wjjThfc| and |w^gfe| in the PR preprocess¬ 
ing stage for the MA stage, and the PR beamforming stage 
for the BC stage, respectively. The two operations provide the 
desired lower bounds for A (i^di, ?x 2 d 2 ) and |w^hfc| (and 
|w^gfc| ), respectively. These desired lower bounds reshape 
the upper bound of the instantaneous end-to-end SER pf E 
in ©, which is optimized relaying on by the RS approach 
of (flOl) to achieve full diversity order. As a result, the actual 
averaged end-to-end performance is also confirmed to exhibit 
full diversity order. 

A. PR strategy 

1) PR Preprocessing in the MA stage: Assume that Rl is 
selected. Considering the DDs A (u\d\, ^ 2 d 2 ), (di,d 2 ) G D , 
which determines of ©, we propose the PR preprocess¬ 
ing to choose u\ and U 2 such that a reasonable lower bound on 
A (i/idi, i/ 2 d 2 ) is given. More specifically, the motivation of 
the PR preprocessing comes from the observation that without 
any preprocessing, when d\ ^ 0 ,d 2 7 ^ 0 , A (i/idi, iz 2 d 2 ) can 
be lower bounded by 

A 2 (u 1 d 1 ,u 2 d 2 ) = ||gfc<fe || 2 + ||h fe di || 2 + 25ft {hf gkd\d 2 u\u 2 } 
= Ilgfcd 2|| 2 + ||h fc L || 2 + 2 \h. k g k d\d 2 \ cos (Zd 2 - Zdi + ip k ) 

> 2 ||gfed 2 || ||h*4i|| - 2 \\g k d 2 \\ \\h k di\\ \cos (Zd 2 - Zdx + <p k )\ 

> 2r kd 2 min (1 - |cos {Zd 2 - Zdi + <p k )\), (11) 

where Lp k = , referred to as the EA of R*, is the angle 

between h fe and g k , T k = min { 11 h fc 11 2 , ||g fe || 2 }, and dmin is 

the minimum value of all possible nonzero values of \di \ and 
|d 2 |. Note that is uniformly distributed over [0,27r]. This 
makes Zd 2 — Zdi + Pk also uniformly distributed. Hence, 
there is a nonnegligible probability that the lower bound in 
CD becomes close to 0. This event then dominates is 
given by ©. 

In order to decrease we eliminate the randomness of 

ipk by judiciously designing u\ and U 2 as 

ui = exp (j (v + pk)) ,i/ 2 = l, (12) 

where j = \/~T> and v is a constant rotation angle that is 
specified later. For the sake of explicit clarity, the geometrical 
interpretation of this PR preprocessing operation is illustrated 
in Fig. 2. Furthermore, the effect of PR is analytically con¬ 
firmed by the following proposition: 
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u l d l h k +u 1 d 2 g k 




Fig. 2. Subfigures (a) and (b) describe the PR preprocessing operation in 
the special case of Zcfi = Zcfe = 0. Subfigure (a) illustrates A (uicfi, 
without any preprocessing in this special case, highlighting how ipk impacts 
A (uidi, U 2 CI 2 ). It is observed from subfigure (b) that once ipk is pre¬ 
canceled, exp (j(pk) h/c<^i and gfccfo coincide along the same direction, hence 
A (u±di , u^d^) is improved. Generally, there always exists (d± , c^) E D 
with Zcfe —Zcfi / 0, thus both <pk and Zcfo—Zcfi are involved in determining 
\(u\di,U 2 d 2 ) as shown in subfigure (c). Subfigure (d) illustrates how 
\(u\di,U 2 d 2 ) is increased by the PR preprocessing in this general case. 
Subfigure (d) shows the rotation of two angles ipk and v, while the rotation 
of pk aims to cope with the angle between and gfc.The rotation of v on the 
other hand intends to adjust the influence of Zcfo — Zd\ on A (uidi, ^ 2 ^ 2 ). 
Finally, comparing subfigures (c) and (d), we see that A (uidi, ^ 2 ^ 2 ) is 
increased by the PR preprocessing operation. 


P roposition 1. With the PR preprocessi ng, we have X lb = 

fcf k |dmin| 2 , i-e., X(uidi,u 2 d 2 ) > fcfk |d mira | 2 for all 
(c?i, 6 ^ 2 ) G D. As a result, p ] jf A is upper bounded as 

pf' 4 < c 4 q ^|c 3 r fe |d min | 2 ^, (13) 

where C 3 - min (1, p min (u)) with p m m (u) = 

min 2 (1 — |cos (Zcfi — Zcfo + u)|), and C 4 = 

(di ,d,2)ED 

Proof: Please see Appendix A. ■ 

As shown in Proposition 1, C 3 = min (1, p m i n (v)) depends 
only on v. Hence, v controls the upper bound of pjf A . 
If Pmin (v) approaches 0 , C 3 = min ( 1 , p m i n (v)) will also 
approach 0 , which may trivialize the upper bound of p^ A 
given in O- Therefore, we constrain v so that p m i n (v) > 0 
to avoid the occurrence of C 3 = 0. Note that v can be easily 
obtained in an off-line manned . 

2) PR Aided Bidirectional Beamforming in the BC 
stage: For the BC stage, we would also like to make 
min | |w^h/c | 2 , |w^g/e| 2 | be determined by T^. This re¬ 
quirement allows min jlw^h/el 2 , Iw^g /^ 2 j and to be 

3 The optimal value of v can be calculated off-line by v = 
arg max p min (V) . 

n'G[0,27r] 


simultaneously optimized by a RS approach regarding T^. 
Correspondingly, both p^ A in © and pf f in ® can be 
improved by the RS approach. According to this requirement, 
we propose a broadcast scheme detailed as follows: 

If Lk = 1, there is no need to design the beamform¬ 
ing vector Wfc. According to for this case, we have 
min ||w^h fc | 2 , |w£g fe | 2 ) = T fc and 

P k,i — CiQ (yC 2 pfiT k ) • 04) 

For the case of Lk ^ 2, the PR beamforming weights can be 
obtained as 


w k 


\ (exp (~j<pk) qf + jqf), (15) 


where ip k = 


h fc s h 


' \h.Pgk | 

Schmidt orthogonalization that is given by 


, qi and q 2 result from the Gram- 


1 . T _ Sfe 


Qi = -h fe , q 2 = 

r 11 




T 22 


(16) 


in which we have ru = ||h fc ||, r 12 = and r 22 = 


&k 


h kS k h 

111 - Ii2 llt 


as follows. 


. The effect of this choice of w k is given 


Proposition 2. For the case of Lk > 2, when w k is chosen 
according to dZ3, we have min || w fch fe | 2 , |wjfg fc | 2 j = \T k . 
Then, a general upper bound of p^f is given by 


p£? < C X Q 


I C 2 pn 


(17) 


Proof: Please see Appendix B. ■ 

Jointly considering (TTtT) and (fl~4l) . we can see that a RS 
approach relying on Tk is capable of lowering the upper bound 
of Pkf • Based on this insight, an RS approach is given as 
follows. 


B. Relay Selection 

Propositions 1 and 2 state that the upper bounds of p^f 
and p^ A are determined by Tk with some modulation-specific 
constants. This insight indicates that using Tk as the selection 
weight throughout both the MA and BC stages reduces p^f 
and p^f A , which in turn improves the overall end-to-end 
performance. Based on this insight, the RS approach in (ITOl) 
becomes 

k = arg max T^. ( 18 ) 

k = l,...,K 

Since Tk = min {11 \ | 2 ,11 gk \ | 2 }, the proposed criterion in 
m constitutes a MaxMin optimization problem. Interestingly, 
we note that the existing MaxMin RS criterion used in single¬ 
antenna relay networks m-m becomes a special case of 
the proposed general MaxMin criterion. 
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IV. Diversity and SER Analysis of the 
PR-MaxMin-RS Scheme 

Achieving full diversity in the end-to-end SER is the main 
design goal of this paper. In particular, the end-to-end di¬ 
versity order is defined as d = lim — } ? g Pe , where Pe 

denotes the average overall end-to-end SER, and it is generally 
limited by the worst among Pm a, Pbc, i, and Pbc, 2 - These 
three constituent error probabilities denote the average SER 
achieved by the relay in the MA stage, and by the sources Si 
and S 2 in the BC stage, respectively. In this section, we first 
elaborate on how the proposed scheme achieves full diversity 
gain. Then, we provides an explicit SER analysis in order to 
further confirm the advantage of the proposed scheme. The 
given analytical results are applicable to the general antenna 
configuration. Additionally, we also discuss two special cases 
which exhibit some some interesting properties. 

A. Diversity Analysis 

In order to present the diversity performance, we first 
investigate the upper bound of the end-to-end SER pf E 
Based on (HTJ and ® , pf 2£ = ± £- =1 (pf A + pf,f) is 
bounded as 

pf 2B < &Q , (19) 

where p is a constant depending on C 2 , C 3 and C 4 . The 
averaged SER is 

Pe = E {pf' E \ < aE [q Ufo FA } < aE {q (MM) } 

k ^ 


where the second inequality of ([2(H) holds owing to the fact 
thatr^^^= max max min {\h k i\ 2 , \gki\ 2 } (see the proof 

k=l,...,Kl = l,...,L k 

in Appendix C) and the property of £ that aE { Q (y//3p £)} oc 

\ (see the proof in Appendix C). The right side of (f20l) that 
implies d = — ^ L. Due to the a priori knowledge that 

d < L, the result d = L is confirmed. 

Remark: We have proved that with the aid of the PR strategy 
is able to achieve full diversity gain. Furthermore, we will 
also show in the sequel that without the PR strategy, the 
full diversity gain cannot be achieved by only applying the 
MaxMin criterion straightforwardly to the general configu¬ 
ration (K ^ 2 , non-binary modulation). The performance 
enhancement obtained in the former scheme comes from the 
fact that the PR preprocessing and PR beamforming operations 
guarantee that the full diversity gain is achieved in the MA 
and BC stages, respectively. 

We also observe that for some special configurations, even 
if the PR preprocessing is not invoked, full diversity gain can 
still be achieved in the MA stage. Hence, the PR preprocessing 
operation and the corresponding overhead spent for feeding 
back u\ can be removed in such scenario. More details can 
be found in the following two propositions. 

Proposition 3. For the scenario with a single multi-antenna 
relay, i.e., K = 1, and arbitrary modulations, there is no need 
to preprocess the transmitted signal of sources for obtaining 


the full diversity gain. In other words, we may have u\ = 
1, U 2 = 1, and the full diversity order d = L can still be 
guaranteed. 

Proof: In this scenario, all of L antennas are gathered in 
one relay. The transmission in the MA stage is equivalent to 
being happened in a 2 x L MIMO system, where it is well 
known that the ML-based MUD in © can obtain full diversity 
without any preprocessing m ■ 

Proposition 4. For the scenario employing BPSK modulation 
and arbitrary antenna configurations, there is no need to pre- 
process the transmitted signal of sources for obtaining the full 
diversity gain. In other words, we may have u\ = 1,^2 = 1, 
the full diversity order d — L can still be ensured. 

Proof: For BPSK modulation, by enumerating all possible 
cases in D, it is readily observed that there is no element of 
{di 0, d 2 / 0} belonging to D. In other words, d\d 2 = 
0 always holds true, which results in A 2 (uidi, 1 x 2 ^ 2 ) = 

I|gfc<^ 2 II 2 + IM 1 II 2 . Then, \(u 1 di,U 2 d 2 ) > \Jv k |d m j „| 2 
holds in spite of any preprocessing. Following the similar 
operation on (l2Qk we see that d = L for BPSK modulation 
even without any preprocessing. ■ 


B. SER Analysis 

In order to further confirm the proposed scheme’s advantage 
in terms of array power gain, we present the SER analysis in 
this subsection. For simplicity, we focus on the scenario where 
all K relays are equipped with the same number of antennas, 
i.e, L\ = L 2 = • • • Lk = £• In order to get Pe , we firstly 
note that Pe < Pm a + \ J2i=i Pbc,i ( 20 ). Additionally, we 
assume that and follow the same distribution, hence 
Pbc ,1 = Pbc, 2 - Therefore, we obtain 

Pe < Pma + Pbc a- ( 21 ) 


It is noted that Pma depends on the particular modulation and 
the corresponding choice of v in the PR preprocessing scheme. 
Hence, considering the MPAM modulation as an example, we 
provide the following proposition about the choice of v and 
the corresponding Pe. 

Proposition 5. For MPAM modulation and the PR prepro¬ 
cessing operation characterized by m, v = ^ is the best 
choice for maximizing A lb- With this particular choice of v, 
Pma and Pbc a are S^ ven as follows: 


Pma < 


E 

(di ,d 2 )eD(s NC ,snc) 


1 

M 2 7T 



IM 1 II 2 

4 sin 2 0 


^71 


M 2 

4 sin 2 0 


dO , 
( 22 ) 


Pbc a ~ 


2(M-1) 

■ So V’ti 

( 3 p 

) de, 

C — 1 

Mtt 

l sin 2 6>(M 2 — l) 

2(M-1) 

/ 0 f V’Tl ( 

f 3 p 


c > 1 

Mtt 

2 sin 2 0(M 2 -1) 


(23) 


where f> 11 (t) is the moment generating function (MGF) of 
71 = /ill hr II , and the definition of f> 11 {f) is f> 11 (t) = 
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^ 1 (t) = Kj2 ( K k 1 ](-l) k (T(t,2k))-2K(K-l)J2 ( K k 2 )(-l) k (2(t,2k+l)-£(t,2k + 2)) 

k =o ' ' k =0 ' ' 


where 


T (t, k) = 


E 


J 0 ,J 2 , ■ --JC -1 
jo + 32 + ■ ■■jc-t = k 


job'2! • • • jc-\\ (£- 1 )! \k + l + tg 


-[ \ ^+SzLo Pz 


2(t,k) = 


E 


(ky.nt-o 1 (h) jl 


3 o, 32 , ■ ■ - 3 c -1 
jo + 32 + . --jc -1 = k 


! \£+Ej = o l oi 


(*+0 


(oc-w 


£ — 1 + Iji 


(£- 1 )! 
(1 + tg) L 


^-i+Efj-o 1 Hi 


s (^) 


1 V +£ (£ + /-!)! 


(24) 


/o°° fn (7) ex P ((7) £7. in which / 7l (7) is t/ze pdf of variable 

71, azrd zY A gzve/t as (E3. 

Proof: Please see Appendix D. ■ 

Summing up (l23t and (l22l) . an upper bound of Pe is 
given by (ETI) . whose accuracy will be further confirmed 
by the numerical results presented later. From the proof of 
Proposition 5, it is noted that the SNR p is always weighted 
by || 11 or || 11 in the instantaneous SER bound. Tak¬ 


ing 


± k\ 

I hr 


for example, we have ||h^ 


= £f=i 


k,l 


where x2 ,..., c independently follow the identical 
distribution. Hence, the fading of C channels characterised 
by i, 2 , • • •, h%. c respectively are averaged by 11 11 . 
This insight can be applied to 11 11 2 similarly. Therefore, an 
improved array power gain can be observed from the SER 
performance achieved because ||h^|| or ||g^|| collects the 
path diversity of C channels. In addition, since the stochastic 

ii 112 

distribution of has been given in Proposition 5, the 

SER result of Proposition 5 may be readily extended to other 
modulations schemes. 

C. Array Power Gain Analysis Based on an Upper Bound of 
Pe 

The above SER analysis implies that the proposed scheme 
enjoys advantages in the achievable array power gain. In this 
subsection, we provide an asymptotic analysis in order to 
understand the achievable array power gain inherent in the 
proposed scheme. Since the stochastic property of the selected 
channel is complex, it is difficult to obtain asymptotic result 
of the actual average SER. Alternatively, we opt for relaxing 
the actual average SER to a desired upper bound of the 
average SER, whose physical significance is very clear. From 
the asymptotic analysis based on the upper bound, we obtain 
further insights into the array power gain and diversity order of 
the proposed scheme. For the sake of simple description of the 
upper bound, we first define T k = Ya=i min {\h k i\ 2 , \g k i\ 2 } 
and k = arg max T k . Then, we have 


r > r~ 

1 k — 1 fc’ 


( 25 ) 


which is obtained by invoking 

T fe = min {11h fe 11 2 , ||g fc || 2 } = min j ^ \h k i\ 2 , E \9kl\ 2 

l 1 = 1 1=1 

-k k 

> y^min{|/i fei | 2 , \gu\ 2 } = f k . (26) 

z=i 

From (l25k the upper bound of the average SER is given by 

Pe = E { P f *} < «E {Q (fpkU) } < aE | Q (ffFfj } . 

(27) 

Then, we focus on the diversity gain and array power gain 
reflected by the upper bound in (l27l) . Using the method in 
m, the upper bound is approximated as 

Pe < aE j Q j « a {G c g)~ Gd , (28) 

where G c characterizes the array power gain, Gd denotes the 
diversity order. In order to obtain the values of G c and Gd, the 
stochastic property of T^ is investigated. Particularly, noticing 
that 

L k 

Pk = max y^min {\h k i\ 2 , \gki\ 2 } , 

k=1 ’~’ K i=i 

Y % p can be treated as the received SNR by applying a two-fold 
diversity technique (detailed later) to a reference lx L SIMO 
P2P system, where the destination’s L antennas are divided 
into K groups and the k -th group consists of Lk antennas. 

The two-fold diversity technique is detailed as follows. 

1) Firstly, for each value of fc, (fc = 1,... K ), Lk diversity 
branches are combined using the maximum ratio com¬ 
bining (MRC) techniques. These branches are mutually 
independent, and the l -th diversity branch is charac- 
terized by min {\h k i\ 2 , \gid\ 2 }- By applying MRC, the 

aggregated channel gain is \JYa=\ min{|/i fci | 2 , \gkf}, 
which is able to provide a diversity order of Gd,k and 
an array power gain of G c ^> 

2) Secondly, the aggregated channel obtained from MRC is 
then treated as a diversity branch, and such K diversity 
















branches are combined using the selection combining 
(SC) technique. These branches are mutually indepen- 
d ent, and the fc-th diversi ty branch is characterized by 
± min {\hki\ 2 ,\gki\ 2 }- Finally, after using the two¬ 
fold diversity technique, the resultant channel gain is 
which provides a diversity order of Gd and an array 
power gain of G c . 

The two-fold diversity technique is termed as the MRC-SC 
scheme in our paper. Following the result of ll22l (Proposition 
4), we have 


K K 

G d = Y J G d ,k = Y, L k = L - (29) 

k= 1 k= 1 

The second equality in (l29l) is based on = L k , which can 
be obtained from classical performance analysis of the MRC 
technique. Again, relying on ll 22 l . G c is expressed as 


G c = 


lK-\ n 


-Z(L + \) 


(30) 


YlkG L c %Z{L k + \) J 

where Z (t) = / 0 °° x l 1 exp (—x) dx , G c ,h is the array 
power gain from combining Lk diversity branches using 
MRC technique. According to the stochastic properties of 
mm{\h k i\ 2 ,\g k i\ 2 } (/ = 1 ,... L k ) and [22), G c ,k is expressed 
as 


Gc,k — 


lLk — 1 


n^Z(l + L k ) 


Z (1 + L k ) 




-l 

L h 


(31) 


Substituting (1311) into ([30 1) . the array power gain G c is ob¬ 
tained. Notice that a{G c fi)~ Gd is an upper bound of Pe- 
Therefore, in terms of achievable SER performance, using the 
proposed scheme in the two-way system is no worse than 
applying the MRC-SC scheme to the reference 1 x L SIMO 
P2P system, where the destination’s L antennas are divided 
into K groups. The antenna partition in the latter reference 
SIMO system is identical to the antenna distribution among 
K relays in the former two-way system. The analytical result 
of (l29l) reveals that the diversity order only depends on the 
total number of antennas. According to (l30t and (1371) . the 
array power gain of the proposed scheme comes from two 
contributing factors: the first one is the MRC-based processing 
detailed in 1 ), and the array power gain offered by this tech¬ 
nique benefits from multiple antennas equipped in each relay. 
The second one is SC-based processing detailed in 2), the SC- 
based processing technique further strengthens the array power 
gain obtained from MRC-processing technique according to 
(l30l) . And the array power gain offered by this technique 
benefits from applying the proposed MaxMin RS scheme to 
K multi-antenna relays. Since the proposed scheme at least 
obtains such two-fold array power gain, the performance of 
the proposed scheme outperforms the other existing schemes, 
which is corroborated by the simulation results provided in 
Section V given later. 


V. Numerical Results And Discussions 

In this section, we compare the SER of the proposed 
PR-MaxMin-RS scheme, the AF based scheme of ifTTIl . the 



■ - - Analytical AF (4,4) 

■ © - Simulated AF (4,4) 

O ' Analytical AF (6,6) 

.Simualted AF (6,6) 

- Analytical PR-MaxMin-RS (4,4), cf.(20) 

-©—Simulated PR-MaxMin-RS (4,4) 

. Analytical PR-MaxMin-RS (6,6), cf.(20) 

-^-Simulated PR-MaxMin-RS (6,6) 


10 12 
SNR 


Fig. 3. The end-to-end SER performance comparison between the proposed 
PR-MaxMin-RS scheme and the AF scheme, when BPSK modulation is 
employed. 


MaxMin-AS scheme of 113 and the O-TR scheme of E2 
by computer simulations. In all simulations, the channel co¬ 
efficients and noises at each antenna are i.i.d. complex-valued 
Gaussian random variables with zero mean and unit variance 
(i.e., cr 2 = 1). In all given figures, the horizontal axis labeled 
by SNR (defined as 5)-) indicates the transmit power of the 
sources. The numerical results of the AF based scheme and 
the O-TR scheme are confirmed by the analytical results given 
in ini and m. Different relay configurations are considered 
as well. 


A. Two Distributed Multi-Antenna Relays 

In Fig. 3 and Fig. 4, the SER performance of the pro¬ 
posed PR-MaxMin-RS scheme, the O-TR scheme, and the 
AF scheme is compared using numerical simulations. Three 
configurations are considered: L\ = L 2 = 4, L\ = L 2 = 6 
and Li = L 2 = 8, which are denoted as (4,4), ( 6 , 6 ), and ( 8 , 8 ), 
respectively. BPSK is employed. The value of p is set to 2 
such that C 2 pmin< \w^h k \ , \w k g k \ > equals to -Af-. As 
shown in both figures, the three schemes considered achieve 
full diversity gain in all configurations. Furthermore, the 
proposed PR-MaxMin-RS scheme attains significantly higher 
array power gains than the AF scheme and the O-TR scheme 
do. This observation can be explained by the fact that the 
proposed PR-MaxMin-RS scheme fully utilizes all antennas in 
the selected relay to combine the signal power throughout the 
MA and BC stages. By contrast, the O-TR and the AF schemes 
only employ one selected antenna among the available L 1 +Z /2 
antennas. This view is also supported by the observation from 
Fig. 3 that the performance gap is enlarged as Li and L 2 
increase. In addition, the analytical upper bound of the SER in 
(I2TT) is also illustrated . It is seen that when the SNR becomes 
high, the upper bound converges to the numerical results. 
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Fig. 4. The end-to-end SER performance comparison between the proposed 
PR-MaxMin-RS scheme and the O-TR scheme, when BPSK modulation is 
employed. 

B. A Single Relay with Multiple Antennas 

In Fig. 5, we consider the centralized scenario where all 
relaying antennas are gathered in a single relay. The MaxMin- 
AS scheme and the transmission beamforming (TB) scheme of 
ED are compared with the proposed PR-MaxMin-RS scheme. 
According to Proposition 3, the PR preprocessing operation 
are removed. Hence, the corresponding overhead is avoided in 
the proposed PR-MaxMin-RS scheme. As shown in Fig.5, the 
TB scheme and the proposed PR-MaxMin-RS scheme achieve 
the similar performance, which is confirmed by the analyt¬ 
ical result as well. Moreover, the numerical and analytical 
SER results of the MaxMin-AS scheme are also illustrated 
in Fig. 5. We can see that the proposed PR-MaxMin-RS 
scheme achieves higher array power gain than the MaxMin- 
AS scheme. This performance gap arises from the fact that the 
MaxMin-AS scheme only exploits the channel associated with 
one antenna in the BC stage, while the proposed PR-MaxMin- 
RS scheme combines the transmitted power of all antennas for 
the sources. 

C. Multiple Distributed Single-Antenna Relays 

In Fig. 6 , we compare the SER performance of the consid¬ 
ered schemes in the configurations where there are 4 and 8 
single-antenna relays. 4PAM modulation is employed and p is 
set to 1. We consider the proposed PR-MaxMin-RS scheme, 
the AF based scheme and the straightforward MaxMin-RS 
scheme dispensing with any PR preprocessing. Two config¬ 
urations, i.e., L = 4 and L = 8 are examined. It is observed 
that the proposed PR-MaxMin-RS scheme achieves the full 
diversity order of L. This is consistent with our analytical 
result given in Section IV. Nevertheless, it is obvious that the 
straightforward MaxMin-RS scheme fails to achieve the full 
diversity. This observation shows that the PR preprocessing 
operation is critical for the MaxMin-RS scheme to achieve 
full diversity. Furthermore, the PR-MaxMin-RS scheme also 
achieves better performance than the AF scheme. This can be 



Fig. 5. The end-to-end SER performance comparison between the proposed 
PR-MaxMin-RS scheme, the MaxMin-AS scheme and the TB scheme, for 
the cases of K = 1, L = 3, 6, when BPSK modulation is employed. 



Fig. 6. The end-to-end SER performance comparison between the proposed 
PR-MaxMin-RS scheme, the straightfoward MaxMin-RS scheme without PR 
preprocessing, and the AF scheme, when K = L = 4, 8, and 4PAM is 
employed. 

explained by the fact that in the AF scheme, the selected relay 
forwards the desired signal together with the undesired noise, 
but the proposed PR-MaxMin-RS scheme only broadcasts the 
intended signal. In addition, we observe that the analytical 
SER given by (f2Th converges to the numerical results in the 
high-SNR region. In summary, although the proposed RS 
scheme optimizes the upper bound of the instantaneous SER, 
as shown in rather than the actual SER, the numerical 
results have shown that the average SER performance exhibits 
full diversity gain as well as better array power gain than 
existing schemes. The numerical results also corroborate our 
analytical results, demonstrating that the upper bound of SER 
given by © is suitable and the tightness of the upper bound 
does not impose a negative impact on the full diversity gain 
and array power gain achieved. 
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Fig. 7. The end-to-end SER performance comparison between the proposed 
PR-MaxMin-RS scheme, the O-TR scheme, and the AF scheme, where K = 
L = 4, and BPSK is employed. 

Fig. 7 provides the SER comparison between the proposed 
PR-MaxMin-RS scheme, the O-TR scheme of (12) and the AF 
scheme of (ED, when BPSK modulation instead of 4PAM is 
employed. This figure shows that the proposed PR-MaxMin- 
RS scheme achieves better performance than the AF scheme, 
but is inferior to the O-TR scheme by no more than ldB. We 
point out that the NCS is assumed to be accurately known 
by relays in the O-TR scheme. By contrast, in the proposed 
PR-MaxMin-RS scheme, the relays do not need to know the 
accurate NCS. Compared with the PR-MaxMin-RS scheme, 
the O-TR scheme assumes that more knowledge about the 
NCS is available to the relays. Therefore, under this condition 
it is reasonable for the O-TR to obtain better performance than 
the PR-MaxMin-RS scheme. Moreover, in the O-TR scheme 
all relays need to decode their received signals in the MA 
stagfl by comparison, such requirement is removed in the 
PR-MaxMin-RS scheme. On the other hand, we acknowledge 
that our proposed scheme requires CSI feedback, while the O- 
TR scheme does not need. Consequently, our proposed scheme 
is not well suitable for the fast-varying fading scenario, where 
the channels only keep static during a small number of symbol 
durations. In principle, the CSI feedback may be employed 
frequently as the channel varies. However, the overhead cost 
for CSI feedback is no longer negligible in this scenario. It 
is important to relax or remove the CSI feedback constraint, 
which will be studied in our future work. 

D. Impact of Channel Estimation Error 

Finally, although our analysis focuses on an RS approach 
under the assumption of perfect channel knowledge, in Fig. 8 
we also show the impact of imperfect CSI encountered during 
the RS process. The channel estimation error is modeled as 
an complex-valued additive white Gaussian noise (AWGN) 

4 Comparing the decoded result with the known NCS, each relay checks 
whether itself should be selected in the BC stage. 



Fig. 8. The end-to-end SER performance of PR-MaxMin-RS scheme under 
different levels of channel estimation error, when K = L = 2, and QPSK is 
employed. 

having a variance S 2 . This AWGN is added to the instan¬ 
taneous channel coefficient and affects the RS performance 
ED- More specifically, Fig. 8 shows the SER performance of 
the proposed PR-MaxMin-RS scheme under different levels 
of channel estimation error indicated by S 2 = 0.1,0.5,0.01, 0. 
The simulation employs QPSK modulation and two single¬ 
antenna relays. It is observed from Fig. 8 that when the 
channel estimation error is serious ( S 2 = 0.1, 0.5), the SER 
performance of our proposed scheme is significantly degraded, 
which also indicates a deteriorated achievable diversity gain. 
As the channel estimation error increases, the SER per¬ 
formance degradation consistently becomes stronger, which 
yields a loss of the achievable diversity gain at high SNRs. 
On the other hand, when the channel estimation error is 
alleviated to S 2 = 0.01, the achievable SER performance and 
the corresponding diversity gain of our proposed scheme are 
close to those of the perfect CSI scenario. This observation 
reveals that our proposed scheme is sensitive to the seri¬ 
ous channel estimation error. Therefore, a high-performance 
channel-estimation algorithm is required to make the proposed 
scheme work effectively. Alternatively, it also implies that 
developing robust RS methods which are insensitive to channel 
estimation error is of great importance in the future work. 

VI. CONCFUSIONS 

This paper investigated a family of PR strategies that 
facilitate the MaxMin-RS criterion to achieve full diversity 
in DF aided TWR systems. Specifically, when the sources 
transmit to the selected relay, they rotate the transmitted 
constellation symbols according to the phases of the selected 
bidirectional channels. When the relay broadcasts the decoded 
NCS, the PR strategy is employed again so that the NCS is 
rotated corresponding to the phases of the broadcast channel. 
By judiciously using the PR strategy in both of the the 
two aforementioned phases, the system eventually achieves 
full diversity in terms of the end-to-end SER. Additionally, 
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we show by both numerical and analytical results that the 
proposed PR-MaxMin-RS schemes outperforms other existing 
relaying schemes under a variety of relay configurations in the 
considered TWR system. 

Appendix A 

The proof of Proposition 1 

Proof: Depending on the value of 

A 2 (uidi, U 2 df) has two kinds of lower bounds. In the 
case of d\ / 0 and 7 ^ 0 , 


where the second equality follows by using = rnqf given 
in (l36l) . and the last equality follows by employing qiq^ = 0 . 
Hence, as shown in ® is given by 

pf.f - C,Q U^\r „ P ) <i> CiQ , 

7 (39) 

where the equality (a) follows by exploiting the fact that 
l|h/c || 2 — r ii as given by (f37l) . Similarly, the received signal 
of S 2 becomes 

VBC, 2 = V^gfe WkSNC + ni 


A 2 (u 1 d 1 ,u 2 d 2 ) = ||gfcM2^i + hkUid 2 \\ 2 

= llgfc^ill 2 + ||h fc d 2 || 2 + 25R{exp(jw)exp(jV fe )gfh fe 4^i} 

= llgfcdill 2 + ||hfcd 2 || 2 + 2 |gf h k d* 2 di\ cos(Zdi - Zd 2 + v) 

> 2 ||gfedi|| ||h fe d 2 || - 2 ||gfedi|| \\h k d 2 \\ |cos(Z<fi - Zd 2 + v)\ 

> r kd 2 min 2 (1 - I cos (Zdi - Zd 2 + u)|), (32) 

'-v-' 

p(v,di,d 2 ) 

where the third equality holds owing to by gjfh/. = 
|g^h/g | exp (— j(pk)' Note that p{v,d\,d 2 ) is a constant de¬ 
pending on v for given d\ and d^ 

On the other hand, for the cases {d\ = 0, d 2 7 ^ 0} and 
{di 7 ^ 0 , d 2 = 0 }, we have 



Or 2 qi + r 22 q 2 ) (exp (-jipk) qf + jqf) s NC + m 
(r 12 exp (-j(fik) + jr 22 ) s NC + rq 


(|t12 | + jV22)SNC +m 


(40) 


where the second equality follows by using g*. = ri 2 qf + 
r 22 q^ as given in(l3Sb, and the last equality follows by <pk = 


h fc sfc 


' l h fgfcl 


= Zri 2 . Then, we have 


Pm = CiQ 



(|ri 2 | + jr 22 ) 


(41) 


A 2 (ixidi,i/ 2 d 2 ) ^ T/cd; 


min’ 


(33) 


Combining (l32b and (l33l) . a general lower bound of 

A (izidi,iz 2 d 2 ) is given by 


= Ci <2 


C 2 p/i 


(|r*i 2 | 2 + 7 * 22 ) = CiQ 


C 2 pp 


Igfcll 


where the second equality is based on the fact that is real- 

= ki 2 | 2 + 


A (uidi, iz 2 d 2 ) ^ min 

(d u d 2 )eD 


valued and the last equality results from | |g^ | 

/ /r W 2 !„(,] ,1 mr tf- \l r 22| 2 as given in <03. Combining d39)i, (SB and ([Hi, we 
\V k mm ' V 1 ’ 2 ’ k mm [obtain the general upper bound for pf? as stated in (fl7l>. 


= JC 3 T k d 2 iri 


(34) 

(35) 


Appendix C 
The proof of (f20li 


where the last equality holds due to C 3 = min ( 1 , p min (O) 
and Pmin M = min p(di,d 2 ,v). Substituting (l34t ) into 

(di ,d 2 )(zD 

©, (fl3b is obtained. Hence, the proof of Proposition 1 is 
established. ■ 


Appendix B 

The proof of Proposition 2 
Proof: According to we also have 

hfc = riiqf, gfc = n 2 qf + r 22 qf. 


and 


|hfc|| 2 = r 2 l5 ||gfc| 


1^12 | 2 + r 


2 

22 * 


(36) 

(37) 


Based on above properties of h/c, g^ and w&, the impact of 
the PR beamforming scheme on is revealed as follows. 
Firstly, Si receives the signal ysc ,l in form of © as 

VBC,1 = V^h^WfcSjVC + Til 

rnqi (exp (~j<pk) qf + jqf) s NC + n x 

rn exp (-j<Pk)sNC + n lf (38) 



Proof: First, let us prove I" k ^ Note that 

T k = max r fe = max min {11h fc 11 2 ,11gfc 11 2 } - (42) 

k=l,...,K k=l,...,K 


For any k E {1,... K}, we have 

H 

r s > r k >J2 min i\ h kf^Ski\ 2 } > min {|/ifcil 2 , Ifliil 2 } = £> 

1=1 

(43) 

where the second inequality follows by invoking 

{ L k L k \ 

i=i i=i J 

L k 

>J2min{\h kl \ 2 ,\g kl \ 2 }. (44) 


1 = 1 


Let 


(XA - argmax min {\h k i\ 2 , \gki\ 2 } . (45) 


Then, (l44b implies that T^ > ( = min {\h ^\ 2 , |g^| 2 }. Next, 
we consider the stochastic property of According to (El), 
the Ith antenna in the kth relay is just the selection result of 
the MaxMin-AS scheme of lITTl applied to all the L relaying 
antennas. It is given in fl3l that aE {Q {y/Pii£) } oc ^ 7 . ■ 

















12 


Appendix D 

The proof of Proposition 5 
Proof: We rewrite (f32l as 


A 2 (uidi,u 2 d 2 ) > T k \d min \ 2 2(l - | cos (Zdi - Zd 2 +t>) |), 

^ ^ ^ 

p(v,di,d 2 ) 

(46) 

where Zdi and Zd 2 equal either 0 or tt when MPAM is 
employed. If v is set to | cos (Zdi — Zd^ + v) | is 0. 
Correspondingly, p (v, di, cfe) always approaches its maximum 
value of 2 . The lower bound of A (uidi, ^ 2 ^ 2 ) in (l32l) is thus 
maximized. Furthermore, by substituting v — | into (l32k the 
average SER in the MA stage becomes 

Pma = E Ak E h .M\/f (iisrfuiMsii 2 ))). 

(di,d 2 )er> k K v v / > 

(47) 


We 


first 


the 




f (l|g^2|| 2 


reshape 

|hfedi || 2 )^ in (l47l) as follows: 


term 


g E h _ |<3 (^f (Ilgrf + IIMxll 2 ))| 


E 



f (llg^2|| 2 


+ llMl|| 2 ) 


2 sin 2 # 



\\di\\ 2 

4 sin 2 0 


^72 




where i/j 7i (t) and ip l2 (t) denote the MGF of 71 = /i||h ^,|| 2 
and 72 = //||g£,|| 2 , respectively. The second equality is 
established on the fact that and are independently 
conditioned on k and that g ^ and h^, have the same distri¬ 
bution, -0 7l (t) = (£). On the other hand, in this case, 

the selected relay has only one antenna. Hence, §nc is 
broadcasted straightforwardly. As a result, Pbc, 1 is given by 


Pbc, 1 = E 



2 (M — 1) /*f 


M 7T 


1: 




f P£mpam\ 
’ 71 V sin 2 0 7 


d<9. 


(49) 


In order to get Pma and Pbc, 1 , we focus on ^ 7l (t) in what 
follows. First, we compute the cdf (Cumulative Distribution 
Function) of ||h ^|| 2 as 


F( 2 ) = Pr{||h s || 2 <4 

= Pr {ll h feH 2 <^|ll h fcH 2 < llgfcl| 2 } Pr {ll h fcH 2 < IlSfcll 2 } 

+ Pr {ll h feH 2 <^|ll h feH 2 > llgfcl| 2 } Pr {ll h fcH 2 > llgfcII 2 } 

= ipr{||h^|| 2 < ^|||h fe || 2 < Hgfell 2 } 

+ ipr{||h fe || 2 ^ ^|||h fe || 2 > llgfcll 2 } - (50) 

The last equality in results from the fact that ||h ^|| 2 
and ||g ^.|| 2 are i.i.d conditioned on k, which implies 
that Pr{||h fc || 2 <||g fc || 2 } = Pr{||h fe || 2 >||g fc || 2 } = \. 
Furthermore, it is noted that Pr{ || m | 2 ^;?||l h fcll 2 ^ ||g fc || 2 } = 
Pr {min {|m|| 2 , ||g fc || 2 } < z), Pr {||h fe || 2 < z\\\h k \\ 2 > ||g ,|| 2 } 
Pr{max{||hj.|| 2 , ||g fc || 2 } < z}. Therefore, d50b becomes 


F ( 2 ) — 2 P min (z) 4” 2 Fmax (z) 5 


(51) 


where Fmin(z) is the cdf of min {||h^|| 2 , ||g^|| 2 }, and 
Fmax (z) is the cdf of max{||m| 2 , ||gfe|| 2 }- From ED, we 
have 


.(*)= 



< £-1 

exp (~z) Y 
K 1=0 



K 


(52) 


On the other hand, Fmax (z) = 

Pr {max{||h fc || 2 , ||g^.|| 2 } < z} . In order to get Fmax (z\ 
we decompose the event {max {||h^,|| 2 , ||g£,|| 2 } ^ zj as 
follows: 


{max{||h^|| 2 , ||gj 2 } < zj 

= U {the/cthrelay is selected, max {IIh/e || 2 , ||gfc 11 2 } 

k=l,...,K L f J J 

= fc=1 u K {uk ^max{||h fe || 2 ,||g fc || 2 },max{||h fc || 2 ,||g fc || 2 } < z) 

(53) 

where u k = max min {||h fe <1| 2 , ||g fc '|| 2 }. Then, 
k'=l,...,K,k'z£k 

Fmax ( 2 ) = Pr {max {||h^.|| 2 , ||g fc || 2 } < 2 } 

K 

= E Pr {“ fc ^ max {11hfe11 2 , ||gfc|| 2 } ,max {||h fc || 2 , ||g fc || 2 } < 2 } 

k=1 
K ~ z 

= E / Pr {uk < ||hfc|| 2 < z,u k < ||gfc|| 2 < z\u k } f (u k ) du k 

k=i Jo 

r z ( 1 1 \ 2 

= K J ( exp (-z) ^2 Jy Z - ex P (~ u ) J2ji u J f ( w ) du ' 

(54) 

where the last equality is based on the fact that ||h ^|| 2 and 
||g ^|| 2 independently follow Gamma distribution conditioned 
on k. Due to the symmetry, the pdf (Probability Distribution 
Function) of u k is identical for all k , we thus remove the 
subscript of / (u k ) to arrive at the last equality. Note that u k 
can be interpreted as the maximum of min{||hfc|| 2 , ||g/c|| 2 } 
among K — 1 values. The pdf of u k can be given by the 
derivative of F m i n (z) while replacing K with K — 1. Based 
on these intermediate results by a few much trivial and tedious 
calculations, the expression of W (t) in (l24l) is obtained. 
With this expression of ^71 (t), we complete the proof of 
Proposition 5. ■ 
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